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1 .  INTRODUCTION 


The  Program  Manager  (PM)  for  BAT  (PM-BAT)  and  his  associated  contractors 
(CAS,  Incorporated,  and  Dynetics,  Incorporated)  have  a  requirement  to  find 
analytical  distributions  to  adequately  represent  the  observed  wind  speed 
distributions  obtained  under  different  meteorological  conditions,  at 
different  time  periods,  and  at  different  geographical  regions  of  the  world. 
These  distributions  are  required  to  test  the  operational  capabilities  of  the 
BAT  system.  The  U.S.  Army  Atmospheric  Sciences  Laboratory  (ASL)  was 
requested  to  perform  the  task. 

Dramatic  changes  in  the  atmosphere  occur  each  day  on  the  average  as  the  sun 
rises  and  sets.  Because  diurnal  variations  are  most  prominent  in  the 
atmosphere ' s  lower  layers ,  changes  in  boundary  layer  winds  have  been 
documented  as  well  as  for  other  atmospheric  variables  (Holton,  1967) . 
Frequency  of  occurrence  of  low  wind  speeds  has  important  applications  in  air 
pollution  analysis  or  for  nuclear-biological-chemical  (NBC)  predictions  for 
military  applications.  The  predictions  of  frequencies  of  high  wind  speeds 
also  have  Important  applications  in  the  areas  of  structural  design  and 
assessment  of  potential  output  from  wind  powered  generators . 

Most  wind  speed  data  used  in  this  analysis  are  from  measurements  made  by 
rotary  cup  anemometers.  For  low  wind  speed  measurements,  two  errors  are 
inherent  in  using  cup  anemometers:  static  friction  and  dynamic  friction. 
A  starting  wind  speed  must  be  sufficient  to  overcome  static  friction  and 
cause  the  cups  to  turn.  Dynamic  friction  becomes  a  factor  as  the  wind 
decreases  from  a  speed  above  the  starting  threshold  until  the  cups  stop 
turning.  Higher  wind  speed  measurements  are  also  subject  to  an  error  known 
as  overspeeding.  Overspeeding  occurs  because  the  anemometer  responds  more 
quickly  to  an  increase  in  speed  than  to  a  decrease  in  speed.  In  a  turbulent 
wind  field  the  mean  value  of  the  anemometer  reading  will  be  higher  than  the 
true  mean.  Also,  an  error  caused  by  vertical  wind  fluctuations  is  of  the 
same  sign  as  the  overspeeding  error.  The  additive  effect  of  these  errors  can 
cause  measurements  to  be  greater  than  the  true  wind  speed  by  as  much  as  10 
percent  (Mage,  1980). 

Investigators  have  examined  various  statistical  models  for  wind  speed 
frequency  distributions  that  range  from  calm  to  the  maximum  observed  wind 
speed.  A  Planck  distribution  was  investigated  by  Dlnkelacker  (1949)  and 
Wentink  (1974).  A  Rayleigh  distribution  (a  chi-square  distribution  with  2 
degrees  of  freedom)  has  been  used  by  Narovlyanskii  (1968)  and  Baynes  (1974) . 
Crutcher  and  Baer  (1962)  showed  that  the  bivariate  normal  distribution  is 
adequate  for  most  wind  samples.  The  Pearson  Type  III  (gamma)  was  used  by 
Putnam  (1948)  and  Sherlock  (1951)  to  describe  wind  speed  distributions.  The 
Weibull  distribution  was  used  to  fit  both  upper  air  data  (Baynes  and 
Davenport,  1975)  and  surface  wind  speed  data  (Wentink,  1976).  Justus  et  al., 
(1976)  were  successful  in  fitting  the  Weibull  distribution  to  data  from  135 
sites.  The  Weibull  model  has  been  an  especially  useful  tool  for  wind  power 
analysis  (Hennessey,  1977).  Stewart  and  Essenwanger  (1978)  examined  the 
Weibull  distribution,  and  Luna  and  Church  (1974)  examined  the  log-normal 
distribution  as  models  for  the  distribution  of  wind  speeds  Independent  of 
direction. 
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2.  CLDUTOLOGY  DATA  BASE  AVAILABLE 


The  ASL  maintains  a  climatology  data  base  containing  historical  weather 
observations  from  more  than  1000  locations  in  North,  Central,  and  South 
America;  Northern,  Central,  and  Southern  Europe;  Mideast  and  North  Africa; 
Southwest,  Southern,  and  Southeast  Asia;  and  Korea.  These  data  were  obtained 
over  a  period  of  years  from  the  U.S.  Air  Force  Environmental  Technical 
Applications  Center  (USAFETAC)  at  Scott  Air  Force  Base,  Illinois.  The 
period- of- record  for  the  data  from  each  station  is  generally  11  years  or  less 
and  may  be  anywhere  from  1966-76  to  1976-86,  depending  upon  which  year  the 
data  were  obtained  from  the  USAFETAC. 

The  regions  of  interest  to  the  PM-BAT  were  the  European  Lowlands,  the 
European  Highlands,  and  the  Hideast  Desert  geographical  areas  (Miers  et  al. , 
1985a;  Miers  et  al.,  1985b).  The  stations  located  in  each  of  these  regions 
are  listed  in  tables  1  through  3.  The  period-of-record  for  the  data  from 
most  of  these  stations  is  1966-76. 

For  the  two  European  regions,  the  PM-BAT  was  interested  in  winter  and  annual 
wind  speed  distributions  obtained  by  grouping  the  observations  from  all  hours 
of  the  day.  Winter  for  these  regions  is  defined  to  be  November,  December, 
January,  and  February.  For  the  Mideast  Desert  region,  the  PM-BAT  was 
interested  in  stammer  and  annual  wind  speed  distributions  obtained  by  grouping 
the  observations  taken  at  night  into  one  group  and  the  observations  taken 
during  the  daytime  into  another  group.  Summer  for  the  Mideast  Desert  region 
is  defined  to  be  June,  July,  and  August. 

The  3 -hourly  observations  (0000,  0300,  ...,  2100  GMT)  in  the  data  base  tend 
to  be  more  reliable  and  consistent  over  time  than  the  observations  taken  at 
the  other  hours  of  the  day.  Therefore,  the  3 -hourly  observations  were  the 
only  observations  included  in  the  analysis  for  any  region,  season,  or  period 
of  the  day. 

3.  WEATHER  CLASSES 

The  PM-BAT  was  interested  in  the  wind  speed  distributions  derived  from 
observations  taken  in  different  weather  conditions  or  classes.  A  weather 
class  is  a  specific  configuration  of  cloud  cover  and  cloud  base  height, 
visibility,  and  obscuration  type.  These  observation  weather  classes  were 
defined  by  representatives  of  the  PM-BA'^'.  The  specific  classes  requested  are 
defined  in  table  4  as  WXCL-1  through  WXCL-6,  WXCL-8,  and  WXCL-9.  The  PM-BAT 
requested  the  wind  speed  information  for  WXCL-1  to  WXCL-6  for  the  European 
regions  and  WXCL-1,  WXCL-8,  and  WXCL-9  for  the  Mideast  region.  The  class 
WXCL-7  was  added  by  the  authors  for  the  European  regions  to  allow  the  PM-BAT 
to  see  the  results  of  fitting  analytical  distributions  to  wind  speeds 
observed  during  fog  conditions  without  regard  to  the  type  of  fog. 
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TABLE  I.  A  LIST  OF  THE  STATIONS  LOCATED  IN  THE  EUROPEAN  LOWLANDS 
GEOGRAPHICAL  REGION 


STATION  STATION  NAME  LATITUDE  LONGITUDE  ELEVATION 

■  SiatBES _ (DEG:MIN)  (DEG:M1N) _ (JBECER&l. 


062400  AMSTERDAM/SCHIPHOL,  NL 

062700  LEEUWARDEN,  NL 

062750  DEELEN,  NL 

062900  TWENTHE/ENSCHED,  NL 

063100  VLISSINGEN,  NL 

063440  ROTTERDAM  AIRPORT,  NL 

063500  GILZE-RIJEN,  NL 

063750  VOLKEL,  NL 

063800  ZUID-LIMBURG/MAAST. ,  NL 

064320  CHIEVRES ,  BX 

064510  BRUSSELS  INTL,  BX 

064560  FLORENNES ,  BX 

064790  KLEIN-BROGEL,  BX 

070150  LILLE/LESQUIN,  FR 

070700  REIMS/CHAMPAGNE,  FR 

071690  ST.  DIZIER/ROBINSON,  FR 

091700  WARNEMUNDE,  DD 

091850  GRIEFSWALD/WIECK,  DD 

093610  MAGDEBURG.  DD 

094880  DRESDEN/KLOTZCHE,  DD 

094920  COTTBUS,  DD 

095540  ERFURT/BINDERSLEBEN,  DD 

100260  HUSUM,  FRG 

100340  EGGEBEK/FLENSBURG,  FRG 

101290  BREMERHAVEN,  FRG 

101470  HAMBURG/FUHLSBUTTEL.  FRG 

102180  AHLHORN,  FRG 

102240  BREMEN,  FRG 

103130  MUNSTER,  FRG 

103140  HOPSTEN,  FRG 

103200  GUTERSLOH,  FRG 

103380  HANNOVER,  FRG 

103840  BERLIN/TEMPLEHOF,  FRG 

104000  DUSSELDORF,  FRG 

104020  WILDENRATH,  FRG 

123300  POZNAN 'LAWICA,  PL 


+52:18 

-004:46 

-0004 

+53:13 

-005:46 

+0001 

+52:04 

-005:53 

+0048 

+52:16 

-006:54 

+0035 

+51:27 

-003:36 

+0010 

+51:57 

-004:26 

-0004 

+51:34 

-004:56 

+0013 

+51:39 

-005:42 

+0022 

+50:55 

-005:46 

+0114 

+50:34 

-003:50 

+0060 

+50:54 

-004:28 

+0055 

+50:14 

-004:39 

+0284 

+51 : 10 

-005:28 

+0060 

+50:34 

-003:06 

+0048 

+49:18 

-004:02 

+0095 

+48:38 

-004:54 

+0139 

+54:11 

-012:05 

+0004 

+54:06 

-013:27 

+0002 

+52:06 

-011:35 

+0079 

+51:08 

-013:46 

+0230 

+51:46 

-014:18 

+0068 

+50:59 

-010:58 

+0315 

+54:31 

-009:09 

+0028 

+54:38 

-009:21 

+0020 

+53:32 

-008:35 

+0007 

+53:38 

-009:59 

+0016 

+52:53 

-008:14 

+0049 

+53:03 

-008:47 

+0003 

+51:58 

-007:36 

+0064 

+52:20 

-007:33 

+0039 

+51:56 

-008:19 

+0070 

+32:28 

-009:42 

+0055 

+52:28 

-013:24 

+0050 

+51:17 

-006:46 

+0045 

+51:07 

-006:13 

+0089 

+52:25 

-016:50 

+0092 
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TABLE  2.  A  LIST  OF  THE  STATIONS  LOCATED  IN  THE  EUROPEAN  HIGHLANDS 
GEOGRAPHICAL  REGION 


STATION 

NUMBER 

LATITUDE 

(DEGrMIN^ 

LONGITUDE 
(DEG: MIN) 

ELEVATION 
^ METERS ) 

065900 

LUXEMBOURG/FINDEL,  BX 

+49:37 

-006:13 

+0491 

071810 

NANCY/OCHEY,  FR 

+48:35 

-005:58 

+0336 

072920 

LUXEUIL/ST.  SAUVEUR.  FR 

+47:47 

-006:21 

+0278 

095460 

KALTENNORDHEIM,  DD 

+50:38 

-010:09 

+0487 

095770 

KARL-MARX-STADT,  DD 

+50:49 

-012:54 

+0357 

104380 

KASSEL,  FRG 

+51:19 

-009:29 

+0158 

105020 

NORVENICH,  FRG 

+50:50 

-006:40 

+0117 

105445 

FULDA  AAF,  FRG 

+50:33 

-009:39 

+0305 

106070 

SPANGDAHLEM,  FRG 

+49:58 

-006:42 

+0365 

106100 

BITBURG,  FRG 

+49:57 

-006:34 

+0375 

106130 

BUCHEL/COCHZM,  FRG 

+50:10 

-007:04 

+0477 

106140 

RAMSTEIN,  FRG 

+49 : 26 

-007:35 

+0237 

106160 

HAHN,  FRG 

+49:57 

-007:15 

+0503 

106590 

KITZINGEN  AAF,  FRG 

+49:45 

-010:12 

+0210 

106870 

GRAFENWOHR,  FRG 

+49:42 

-011:57 

+0414 

107080 

SAARBRUCKEN/ENSHEIM,  FRG 

+49:13 

-007:07 

+0322 

107140 

ZWEIBRUCKEN,  FRG 

+49:13 

-007:25 

+0343 

107380 

STUTTGART/ECHTER,  FRG 

+48:41 

-009:12 

+0396 

107630 

NURNBERG,  FRG 

+49 : 30 

-011:05 

+0319 

108570 

LANDS BERG/LECH,  FRG 

+48:04 

-010:54 

+0623 

108580 

FURSTENFELDBRUCK,  FRG 

+48:12 

-011:16 

+0518 

108660 

MUNCHEN/RIEM,  FRG 

+48:08 

-011:43 

+0529 

110100 

LINZ/HORSCHING,  OS 

+48:14 

-014:11 

+0297 

111500 

SALZBURG,  OS 

+47:48 

-013:00 

+0430 

114060 

CHEB,  CZ 

+50:05 

-012:24 

+0474 

115180 

PRAGUE/RUZYNE,  CZ 

+50:06 

-014:17 

+0380 

125100 

SNEIZKA,  PL 

+50:44 

-015:44 

+1603 
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TABLE  3.  A  LIST  OF  THE  STATIONS  LOCATED  IN  THE  MIDEAST  DESERT 
GEOGRAPHICAL  REGION 


STATION 

NUMBER  STATION  NAME _ 

171700  VAN,  TU 

171950  KAYSERI/ERKILET,  TU 

172400  ISPARIA,  TU 

172800  DIYARBAKIR,  TU 

400070  ALEPPO/NEIRAB,  SY 

400160  HASSAKAH.  SY 

400390  RAQQA,  SY 

400610  PAIMYRA/TUDMUR.  SY 

400800  DAMASCUS  NEW  INTL,  SY 

402600  HOTEL  FIVE,  TJ 

402700  AMMAN/KING  ABDULLAH,  TJ 

403100  MA'AN,  TJ 

403560  TURAIF/AL  rJRAYF,  SD 

403620  RAFHA,  SD 

403750  TABOUK,  SD 

403940  HAIL,  SD 

404300  MEDINA/Ai.  MADINAH,  SD 

404380  RIYADH  INTL,  SD 

406080  MOSUL/MOSSOUL,  IQ 

406340  HADITHA/HAQLANIYA,  IQ 

406420  RUTBAH,  IQ 

406500  BAG*.DAD  INT  (ORBW) ,  IQ 

406740  SEMAWA,  IQ 

407060  TABRIZ,  IR 

407430  SABZEVAR,  IR 

407470  SANANDAJ ,  IR 

407540  TEHRAN/MEHRABAD,  IR 

407910  TABAS,  IR 

408000  ESFAHAN,  IR 

408290  ZABUL,  IR 

408410  KERMAN,  IR 

408590  FASA,  IR 

408790  IRANSHAHR,  IR 

409110  MAZARI- SHARIF,  AH 

409130  KUNDUZ,  AH 

409380  HERAT,  AH 

409740  FARAH,  AH 

409860  ZARANJ/NIMROZE,  AH 

409900  KANDAHAR,  AH 

417100  NOKKUNDI,  PK 

417120  DALBANDIN,  PK 

417390  PANJGUR,  PK 

623060  MERSA  MATRUH,  UB 

623180  ALEXANDRIA/NOUZHA,  UB 


LATITUDE 

rDEG:MIN> 

LONGITUDE 

rDEG:MIN) 

ELEVATION 
f METERS) 

+38:27 

-043:19 

+1667 

+38:47 

-035:29 

+1052 

+37:45 

-030:33 

+0997 

+37:54 

-040:12 

+0686 

+36:11 

-037:13 

+0390 

+36 : 30 

-040:45 

+0295 

+35:56 

-039:01 

+0245 

+34:33 

-038:18 

+0395 

+33:25 

-036:31 

+0610 

+32:12 

-037:08 

+0668 

+31:59 

-035:59 

+0767 

+30:10 

-035:47 

+1069 

+31:41 

-038:40 

+0827 

+29:38 

-043:29 

+0443 

+28:22 

-036:35 

+0771 

+27:31 

-041:44 

+0992 

+24:33 

-039:43 

+0646 

+24:42 

-046:44 

+0624 

+36:19 

-043:09 

+0223 

+34:04 

-042:22 

+0140 

+33:02 

-040:17 

+0615 

+33:14 

-044:14 

+0034 

+31:18 

-045:16 

+0006 

+38:08 

-046:15 

+1367 

+36:13 

-057:40 

+0941 

+35:20 

-047:00 

+1538 

+35:41 

-051:19 

+1204 

+33:36 

-056:54 

+0691 

+32:37 

-051:40 

+1598 

+31:01 

-061:29 

+0487 

+30:15 

-056:58 

+1748 

+28:58 

-053:41 

+1382 

+27:12 

-060:42 

+0566 

+36:42 

-067:12 

+0378 

+36:40 

-068:55 

+0433 

+34:13 

-062:13 

+0964 

+32:22 

-062:11 

+0700 

+31:06 

-061:56 

+0478 

+31:30 

-065:51 

+1010 

+28:49 

-062:45 

+0679 

+28:53 

-064:24 

+0848 

+26:58 

-064:06 

+0968 

+31:20 

-027:13 

+0028 

+31:12 

-029:57 

-0003 
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Table  3 .  (cent) 


STATION 

NUMBER 

LATITUDE 

(DEG:MIN) 

LONGITUDE 
(DEG: MIN) 

ELEVATION 

(METERS) 

623330 

PORT  SAID/EL  GAMIL,  UB 

+31:17 

-032:14 

+0003 

623660 

CAIRO  INTL,  UB 

+30:08 

-031:24 

+0112 

6238/0 

EL  MINYA/MINIA,  UB 

+28:05 

-030:44 

+0039 

624140 

ASWAN/ASWAN  DAM,  UB 

+23:58 

-032:47 

+0200 

624230 

FARAFRA,  UB 

+27:03 

-027:58 

+0090 

626500 

DONGOLA,  SU 

+19:10 

-030:29 

+0226 

627210 

KHARTOUM,  SU 

+15:36 

-032:3^ 

+0382 

627300 

KASSALA,  SU 

+15:28 

-036:24 

+0500 

627810 

EN  NAHUD,  SU 

+12:42 

-028:26 

+0565 

TABLE 

4.  DEFINITIONS  OF  THE  NINE  OBSERVATION  HEATHER  CLASSES  IN 

TERMS  OF  CLOUD  COVER.  CLOUD  BASE  HEIGHT.  OBSCURATION 

TYPE,  AND  VISIBILITY 

LABEL 

WXCL-1 

No  Clouds  (no  clouds  or  no  clouds  with  bases  <-  3000  m) 

WXCL-2 

Clouds  with  bases  <-  3000  m 

WXCL-3 

Clouds  with  bases  <-  3000  m  and  Precipitation 

WXCL-4 

Clouds  with  bases  <-  3000  m  and  Advective  Fog 

WXCL-5 

Clouds  with  bases  <-  3000  m  and  Precipitation  and 

Advective  Fog 

WXCL-6 

Radiative  Fog 

WXCL-7 

Fog 

WXCL-8 

Sandstorms  and  Dust 

WXCL-9 

Low  Visibility  (3  <-  visibility  <-  10  km) 

By  World  Meteorological  Organization  (WHO)  standards,  a  fog  can  be  reported 
only  when  the  visibility  is  not  more  than  1000  m.  The  PM-BAT  wished  to  make 
a  distinction  between  radiative  and  advective  fog  conditions.  Observations 
in  the  ASL  climatology  data  base  do  contain  fog  reports,  but  they  do  not 
contain  information  about  the  type  of  fog  reported.  To  satisfy  the  require¬ 
ments  of  the  PM-BAT,  definitions  of  radiative  and  advective  fogs  were  taken 
from  the  USAFETAC  electro-optical  (EO)  climatology  software  (USAFETAC,  1991). 
A  reported  fog  was  classified  as  radiative  fog  if  the  wind  speed  was  less 
than  3.5  m/s  (7  kn) .  Similarly,  a  reported  fog  was  classified  as  advective 
fog  if  the  wind  speed  was  at  least  3.5  m/s. 
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4.  ANALYTICAL  DISTRIBUTIONS 


Three  analytical  distributions  were  selected  for  potential  use  in  repre¬ 
senting  the  observed  wind  speed  distributions.  They  were  chosen  based  upon 
the  work  of  previous  researchers  referred  to  in  section  1  and  the  request 
from  the  PH-BAT  to  specifically  Investigate  the  gamma  distribution.  These 
three  candidate  distributions  were  the  gamma,  log-normal,  and  Veibull.  Each 
of  these  distributions  performed  reasonably  well  at  representing  the  observed 
wind  speed  distributions  for  most  of  the  weather  classes. 

The  gamma  probability  distribution  is  defined  by 


_ _ _ E _ 

pr(a) 


(1) 


where  r(a)  is  the  gamma  function  with  argtiment  a,  x  >  0,  a  >  0,  and  /J  >  0. 


The  log-normal  probability  distribution  is  defined  by 


_ 2 _ ff _ 

xo>/2w 


(2) 


where  n  is  the  expected  value  and  a  is  the  standard  deviation  of  ln(x)  and 

X  >  0. 


The  Weibull  probability  distribution  is  defined  by 


f(x;a,P)  =  (  «)  (^)«-iexp[-(^)«] 
P  P  P 


(3) 


where  x  >  0,  a  >  0,  and  P  >  0. 

In  this  report  n  and  p  are  called  scale  parameters  because  changing  their 
values  either  tends  to  move  the  center  or  skew  the  distribution  toward  hl^er 
or  lower  values .  a  and  a  are  called  shape  par^uaeters  because  changing  their 
values  tends  to  change  the  shape  of  the  distribution  to  either  a  more  peaked 
or  a  flatter  form. 

Three  procedures  for  computing  values  for  a  and  P  or  n  and  a  were  considered 
for  fitting  the  analytical  distributions  to  the  observed  wind  speed  distri¬ 
butions  (histograms).  These  procedures  are  the  method  of  moments,  swximum 
likelihood,  and  nonlinear  least  squares  fit  to  the  observed  cumulative 
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probability  distribution.  For  each  observed  wind  speed  distribution,  the  PM- 
BAT  was  interested  in  determining  the  best  of  the  three  analytical  distri¬ 
butions  and  the  associated  parameter  values.  Once  the  name  of  the  best 
analytical  distribution  and  the  distribution  parameter  values  are  known,  the 
estimation  of  the  tail  probability  of  the  wind  speed  is  possible. 

The  "goodness-of-fit"  of  the  analytical  distributions  was  compared  by 
computing  both  the  mean  absolute  difference  and  the  root-mean-square  (rms) 
difference  between  the  observed  and  each  of  the  three  analytical  cximulatlve 
distributions.  After  thoroughly  comparing  the  analytical  distributions  using 
both  criteria,  we  decided  to  use  the  minimum  rms  difference  to  make  the  final 
selection. 

When  this  criterion  is  used,  the  method  of  moments  generally  yields  parameter 
estimates  that  do  not  perform  as  well  as  the  other  two  procedures,  me  maxi¬ 
mum  likelihood  procedure  tends  to  yield  better  estimates  of  the  parameters 
than  the  method  of  moments.  The  nonlinear  least  squares  fit  gives  the  best 
parameter  estimates,  as  expected  by  the  selection  criterion. 

4 , 1  Method  of  Moments 

The  expected  value  of  x  for  the  gamma  distribution  is  given  by 


E[x]  *  aP  W 

and  the  variance  of  x  is  given  by 

VARlx]  »=  ap*  .  (5) 

These  two  equations  can  be  used  to  compute  values  of  a  and  p  for  the  gamma 
distribution  using  observed  data  to  estimate  E[x)  by  the  mean  value  of  the 
observed  data  and  VAR[x]  by  the  sample  variance  of  the  observed  data  giving 

P  =  (6) 

i?tx3 

and 

a  =  ,  (7) 

This  procedure  for  obtaining  values  for  the  distribution  parameters  is  called 
the  method  of  moments. 

Using  the  method  of  moments  to  compute  values  for  n  and  a  in  the  log-normal 
distribution  gives 

^  =  E[ln{x)  ] 
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and 


a  =  v'VAJtllnUJJ 


(9) 


where  E[ln(x)]  and  VAR[ln(x)]  are  estimated  from  the  data  by  using  the  mean 
value  and  sample  variance  of  the  natural  logarithm  of  the  wind  speed  data. 

The  advantage  of  using  E[ln(x)]  and  VAR[ln(x)]  rather  than  E[x]  and  VAR[x] 
to  estimate  ft  and  a  is  quite  obvious  when  one  considers  the  relationship 
between  the  distribution  parameters  and  E[x}  and  VAR[x]  given  by 

E[x]  =  exp[^  + 

A 


and 


VAR{x\  =  exp (2^  +  o*)  lexp(o^)  -  1]  . 

The  expected  value  of  x  for  the  Weibull  distribution  is  given  by 

E[x]  =  lr(A) 

ft  ft 


(11) 


(12) 


and  the  variance  of  x  is  given  by 


VAR[x]  =  S1[2T(^)  -A[r(i:)]2].  (13) 

CL  CL  CL  CL 

It  is  rather  difficult  to  use  the  method  of  moments  to  compute  values  of  a 
and  p  for  the  Weibull  distribution  using  estimates  of  E[x]  and  VAR[x] 
obtained  from  observed  data.  Considering  the  computational  time  required  to 
iteratively  solve  this  system  of  equations  and  the  general  inferior  quality 
of  the  parameter  estimates,  this  procedure  was  not  used  to  compute  the 
Weibull  distribution  parameters. 

4.2  Maximum  Likelihood 

The  maximum  likelihood  procedure  involves  finding  the  parameter  values  that 
maximize  the  likelihood  function  of  the  distribution.  For  the  gamma 
distribution  these  parameter  values  are  obtained  by  iteratively  solving  the 
following  system  of  equations: 


ln(p)  +  Y(a)  =  E[ln(x)  ] 


(14) 


ap  =  E[x] 


(15) 


where  ♦(a)  is  the  digamma  function,  ln[r(a)]'  ('  indicating  first  derivative 
with  respect  to  a),  with  r(Q)  being  the  gamma  function  with  argument  o. 
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The  maximuin  likelihood  estimates  of  the  log-normal  distribution  parameters 
are  the  same  as  those  given  in  equations  (8)  and  (9)  by  the  method  of 
moments . 


The  maximum  likelihood  estimates  of  the  Weibull  distribution  parameters  are 
obtained  by  iteratively  solving  the  following  system  of  equations. 


g[x*ln(x)  ]  _  _1 
E[x^]  a 


E[ln{x)] 


(16) 


p  -  B[x«]  •  . 


(17) 


4.3  Nonlinear  Least  Squares  Fit 

The  least  squares  fit  was  performed  on  observed  cumulative  distribution 

values  {Pj,  i-1,  2 . N)  corresponding  to  discrete  values  of  the  wind 

speed,  X,  denoted  by  {X^,  i-1,  2,  ...  ,  N),  where  N  is  the  number  of  wind 
speed  intervals.  Each  wind  speed  interval  had  a  width  of  0.5  m/s.  N  was 
determined  by  choosing  the  first  interval  of  the  observed  wind  speed 
frequency  of  occurrence  histogram  such  that  all  intervals  greater  than  the 
Nth  contained  no  occurrences  of  wind  speed. 

The  three  analytical  distributions  are  clearly  nonlinear  with  respect  to  the 
distribution  parameters  a,  fi,  and  a.  To  determine  the  least  squares 
estimates  of  these  parameter  values,  a  procedure  was  used  to  linearize  the 
distributions  (at  least  in  the  neighborhood  of  any  specific  values  of  the 
parameters) . 

Consider  the  integral  of  the  probability  distribution  given  by 

X 

F(2f;a,P)  =  Jf  (x;a,  p)  dx 

0 

By  letting 

a  =  ttp  +  6a 

and 

P  =  Po  +  6P  ,  (20) 


(18) 

(19) 


equation  (18)  becomes 


X 

F(X:aQ  +  6a,  Pg  +  6P)  =  Jf(x;ag  +  6a,  Pg  +  dfi)  dx  .  (21) 

0 

To  a  first-order  approximation  using  small  values  for  6a  and  6/3 
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where 


m!«,p)  -  FU;«„Po>  ‘ 


dF{X;  ttp,  Pq)  _  f  df(x;a,  P) 
da  J  da 


«  “  «o 


aF(X; ao,Po)  _  fdf{x}a,fi) 

ap  "  J  ap 


fi  -  Po 


Equation  (22)  is  linear  In  6a  and  6fi. 


Let  Q  represent  the  sum  of  the  squares  of  the  differences  between  the 
analytical,  FCX^;  a,  /9) ,  and  observed  cumulative,  P^,  distributions  given  by 


Replacing  FCX^;  a,  /3)  by  the  expression  in  equation  (22)  gives 


Taking  the  partial  derivatives  of  Q  with  respect  to  3a  and  6fi  and  setting 
each  derivative  equal  to  zero  to  find  a  minimum  of  Q  give 

E  *  E  =  E 
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and 


1 1 ^  t  -tiPi-  FoU,)]i£^ 


ap 


ap 

(28) 


where  the  terms  involving  Fo(Xi)  and  its  partial  derivatives  denote  the 
function  F(x;  a,  j9)  and  its  partial  derivatives  evaluated  at  X^,  Oq,  and 
These  two  equations  are  linear  in  5a  and  6^. 


The  nonlinear  least  squares  fit  procedure  begins  by  getting  initial  values 
for  Oq  and  /Sq;  solving  equations  (27)  and  (28)  for  6a  and  Sfi  using  the 
definitions  in  equations  (18),  (23),  and  (24);  and  getting  updated  values  of 
a  and  ^  according  to  the  Iterative  equations 


and 


®JC  = 

(29) 

Pjc  ~  Pjc-1  sp 

(30) 

for  k  -  1 ,  2 ,  ...  until  6a  and  6/9  become  very  small  and  convergence  is 
obtained.  The  initial  values  of  the  parameters  may  be  selected  by  any 
method.  Reasonable  initial  values  are  obtained  by  using  the  method  of 
moments  to  get  a^  and  Pq. 

5.  WIND  SPEED  DISTRIBUTIONS 

All  three  of  the  analytical  distributions  considered  in  this  analysis  have 
zero  probability  of  occurrence  at  x  -  0.  The  observed  wind  speed,  however, 
may  be  calm  or  zero.  To  account  for  this  phenomenon  in  all  weather  classes 
except  WXCL-4  and  WXCL-5  (advective  fog  cases),  the  wind  speed  cumulative 
distribution  is  defined  by  the  equation 


W{X;p^,a.P)  =PoS(0)  +  (1  -  Po)F(2f;a,  P)  (31) 

where  X  is  the  wind  speed,  6(0)  is  the  Dirac  delta  function,  Po  is  the 
probability  of  calm  wind  speeds,  and  F(X;  a,  is  the  appropriate  analytical 
cumulative  distribution  function. 

For  weather  classes  WXCL-4  and  WXCL-5,  the  wind  speed  cumulative  distribution 
is  defined  by  a  modification  of  the  previous  equation  to  give 


W(X;Po,a,fi)  =  (1  -  Po)F(X  -  3;a,p)  (32) 

where  the  wind  speed,  X,  is  at  least  3.5  ra/s.  The  probabilities  of  calm  wind 
speeds  for  each  region,  time  period,  and  weather  class  combination  are  listed 
in  table  5.  Similarly,  the  numbers  of  noncalm  wind  speed  observations  are 
listed  in  table  6. 
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TABLE  5.  FROBABIUTIES  OF  CALM  WIND  SPEEDS  FOR  ALL  WIND  SPEED 
DISTRIBUTION  CASES 


Nideast  Desert 


Anniial  Day 


WXCL-1 

WXCL-8 

WXCL-9 


17.75504 

6.83641 

16.04527 


TABLE  6 .  THE  NUMBER  OF  NONCALM  HIND  SPEED  OBSERVATIONS  FOR  ALL 
HIND  SPEED  DISTRIBUTION  CASES 


Regjop _ Season _ Time  Period  Wx  Class  Noncalm  Count 


European  Lowlands 

Winter 

All  Hours 

WXCL-1 

54211 

WXCL-2 

230939 

WXCL-3 

73280 

WXCL-4 

3621 

WXCL-6 

10775 

WXCL-7 

14528 

European  Lowlands 

Annual 

All  Hours 

WXCL-1 

222626 

WXCL-2 

642539 

WXCL-3 

171815 

WXCL-4 

5456 

WXCL-6 

21990 

WXCL-7 

27736 

European  Highlands 

Winter 

All  Hours 

WXCL-1 

31838 

WXCL-2 

182433 

WXCL-3 

62402 

WXCL-4 

4753 

WXCL-5 

217 

WXCL-6 

9705 

WXCL-7 

14502 

European  Highlands 

Annual 

All  Hours 

WXCL-1 

137153 

WXCL-2 

504420 

WXCL-3 

141439 

WXCL-4 

9842 

WXCL-5 

277 

WXCL-6 

16089 

WXCL-7 

26019 

Mideast  Desert 

Summer 

Night 

WXCL-1 

44780 

WXCL-8 

1451 

WXCL-9 

18979 

Mideast  Desert 

Summer 

Day 

WXCL-1 

84355 

WXCL-8 

8081 

WXCL-9 

39687 

Mideast  Desert 

Annual 

Night 

WXCL-1 

138562 

WXCL-8 

6534 

WXCL-9 

84275 

Mideast  Desert 

Annual 

Day 

WXCL-1 

193413 

WXCL-8 

24516 

WXCL-9 

133441 
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The  observed  probability  distribution  was  obtained  by  fonsing  histograns  of 
frequency  of  occurrence  of  observed  wind  speed  with  all  histogram  interval 
widths  being  O.S  m/s.  The  histogram  frequencies  of  occurrence  were  converted 
to  probabilities  and  assigned  to  the  wind  speeds  at  the  midpoints  of  the 
intervals.  The  corresponding  analytical  probability  distribution  was 
obtained  by  evaluating  the  analytical  cumulative  probability  distribution  at 
the  wind  speeds  corresponding  to  the  beginning  and  end  of  each  histogram 
Interval  and  using  the  difference  between  the  two  values  assigned  to  the  wind 
speed  at  the  midpoint  of  the  interval. 

The  mean  absolute  deviation  between  the  observed  and  analytical  cumulative 
distributions  for  each  region,  time  period,  and  weather  class  combination  is 
listed  in  table  7.  Similarly,  the  rms  deviation  between  the  cumulative 
distributions  is  listed  in  table  8.  Generally,  both  the  mean  absolute 
deviations  and  the  rms  deviations  are  small,  indicating  the  analytical 
distributions  adequately  represent  the  observed  wind  speed  distributions. 
The  names  of  the  corresponding  best  distribution  are  listed  in  table  9. 

6 .  ANALYSIS  RESULTS 

6 . 1  European  Lowlands 

The  observed  noncalm  wind  speed  and  analytical  distributions  along  with  the 
observed  and  analytical  cumulative  distributions  for  the  European  Lowlands 
are  depicted  in  figures  1  through  12.  The  figures  correspond  to  winter  and 
annual  time  periods  for  weather  classes  WXCL-1  to  WXCL<4,  VKCL>6,  and  WXCL-7. 
Weather  class  WXCL-5  had  an  insufficient  number  of  observations  for  analysis. 
The  odd  niimbered  figures  depict  the  observed  and  analytical  distributions, 
and  the  even  number  figures  depict  the  observed  and  analytical  cumulative 
distributions.  In  each  figure  the  differences  between  the  observed  and 
analytical  values  are  also  plotted.  The  observed  distributions  are  plotted 
using  a  dashed  line,  analytical  distributions  with  a  solid  line,  and  the 
difference  between  them  with  a  dotted  line. 

The  figures  depict  only  the  noncalm  part  of  the  wind  speed  distributions. 
The  number  of  noncalm  wind  observations  used  to  fit  the  analytical 
distributions  to  the  observed  wind  speed  distributions  is  displayed  in  each 
figure.  In  addition,  the  best  analytical  probability  distribution  and  the 
corresponding  scale  and  shape  parameters  are  listed  in  each  figure. 

Of  the  12  different  weather  class  and  season  combinations  there  were  3  cases 
where  the  gamma  distribution  was  the  best,  3  for  the  log-normal  and  6  for  the 
Weibull  (reference  table  9).  The  log-normal  distribution  was  the  best  for 
weather  class  WXCL-4  (advective  fog  case)  and  the  annual  wind  speed 
distribution  for  WXCL-7  (fog).  The  gamma  distribution  was  the  best  for  the 
winter  wind  speed  distributions  for  WXCL-1,  WXCL-2,  and  WXCL-7.  Overall,  the 
analytical  distributions  perform  well  at  describing  the  observed  wind  speed 
distributions.  The  worst  performance  was  for  the  weather  classes  involving 
fog.  Even  for  those  cases  the  analytical  distributions  perform  fairly  well. 
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6.2  European  Highlands 

Figures  13  through  26  depict  the  same  type  of  wind  speed  Information  for  the 
European  Highlands.  For  this  region  weather  class  WXCL-5  had  enough  noncalm 
wind  speed  observations  to  be  analyzed.  Of  the  14  different  weather  class 
and  season  combinations  there  were  2  cases  where  the  gamma  distribution  was 
the  best,  10  for  the  log-normal  and  2  for  the  Weibull  (reference  table  9). 
The  gamma  distribution  was  best  for  the  winter  wind  speed  distributions  for 
WXCL-5  (advective  fog  and  precipitation)  and  WXCL-6  (radiative  fog).  The 
Weibull  distribution  was  the  best  for  weather  class  WXCL-4  (advective  fog) . 
Overall,  the  analytical  distributions  performed  fairly  well  for  all  cases 
except  WXCL-4  and  WXCL-7.  Again,  the  analytical  distributions  for  weather 
classes  involving  fog  proved  to  have  the  worst  performance,  with  WXCL-4  and 
WXCL-7  being  unsuited  for  analytical  distribution  representation. 

6 . 3  Mideast  Desert 

Figures  27  through  38  depict  the  same  type  of  wind  speed  information  for  the 
Mideast  Desert  region,  except  that  the  distributions  are  for  weather  classes 
WXCL-1,  WXCL-8,  and  WXCL-9  and  the  time  periods  are  summer  and  annual  for 
night  and  daytime.  Of  the  12  different  weather  class  and  time  period 
combinations  there  were  2  cases  where  the  gamma  distribution  was  the  best  and 
10  for  the  Weibull  (reference  table  9).  The  gamma  distribution  was  the  best 
during  the  annual  night  period  for  WXCL-1  (no  clouds)  and  WXCL-9  (low 
visibility) .  For  the  Mideast  the  observed  wind  speed  distributions  are  more 
"lumpy"  than  those  for  the  European  regions.  However,  the  analytical 
distributions  appear  to  adequately  represent  the  observed  distributions.  The 
worst  representations  occur  for  the  daytime  dust  storm  conditions  (WXCL-8) 
when  the  wind  speeds  tend  to  be  greater  and  more  variable. 

7 .  CONCLUSIONS 

The  observed  wind  speed  distributions  are  adequately  represented  by  the 
analytical  distributions  except  in  a  few  cases.  For  these  exceptions, 
weather  classes  WXCL-4  and  WXCL-7  for  the  European  Highlands,  the  observed 
wind  speed  distributions  do  not  lend  themselves  to  analytical  representation 
with  enough  accuracy  to  be  useful.  The  PM-BAT  Intends  to  use  the  analytical 
distributions  to  generate  wind  speeds  by  pseudorandom  numbers.  For  the  two 
exception  cases ,  the  observed  wind  speed  empirical  distributions  should  be 
used  rather  than  the  analytical  distributions  to  generate  random  wind  speeds . 
In  all  cases  whether  analyti'  .1  or  empirical  distributions  are  used, 
generated  wind  speeds  more  than  40  m/s  should  be  discarded.  Furthermore, 
occurrences  of  such  large  values  of  wind  speed  may  indicate  a  problem  with 
the  pseudorandom  number  generator. 
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TABLE  7.  THE  MEAN  ABSOLUTE  DEVIATION  BETWEEN  THE  OBSERVED  AMD 
ANALYTICAL  CUMULATIVE  DISTRIBUTIONS 


Rygivn _ Ss.ft?9D _ T1m«  Period  Vx  Cla«« _ M<«n  Ab»  Dev  (%'i 


European  Lowlands 

Winter 

All  Hours 

WXCL-1 

0.29212 

WXCL-2 

0.29047 

WXCL-3 

0.30860 

WXCL-4 

0.38378 

WXCL-6 

3.47584 

WXCL-7 

0.43066 

European  Lowlands 

Annual 

All  Hours 

WXCL-1 

0.36995 

WXCL-2 

0.39851 

WXCL-3 

0.32126 

WXCL-4 

0.31072 

WXCL-6 

3.28172 

WXCL-7 

0.32688 

European  Highlands 

Winter 

All  Hours 

WXCL-1 

0.33095 

WXCL-2 

0.55635 

WXCL-3 

0.77190 

WXCL-4 

3.59613 

WXCL-5 

1.78300 

WXCL-6 

3.00902 

WXCL-7 

2.29373 

European  Highlands 

Annual 

All  Hours 

WXCL-1 

0.29323 

WXCL-2 

0.49851 

WXCL-3 

0.65162 

WXCL-4 

2.82679 

WXCL-5 

1.80694 

WXCL-6 

2.81105 

WXCL-7 

2.29231 

Hldeast  Desert 

Summer 

Night 

WXCL-1 

0.46519 

WXCL-8 

0.58164 

WXCL-9 

0.46327 

Mideast  Desert 

Summer 

Day 

WXCL-1 

0.44986 

WXCL-8 

0.97730 

WXCL-9 

0.38925 

Mideast  Desert 

Annual 

Night 

WXCL-1 

0.38754 

WXCL-8 

0.47509 

WXCL-9 

0.36817 

Mideast  Desert 

Annual 

Day 

WXCL-1 

0.47917 

WXCL-8 

0.73087 

WXCL-9 

0.40136 
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TABLE  8.  THE  ROOT -MEAN -SQUARE  DEVIATION  BETWEEN  THE  OBSERVED  AND 
ANALYTICAL  CUMULATIVE  DISTRIBUTIONS 


- SMAfin - limft  Period  Wx  Cl>««  RMS  Dev  (%') 


European  Lowlands 


European  Lowlands 


European  Highlands 


European  Highlands 


Mideast  Desert 


Mideast  Desert 


Mideast  Desert 


Mideast  Desert 


Winter 

All  Hours 

WXCL-1 

0.58660 

WXCL-2 

0.54364 

WXCL-3 

0.51059 

WXCL-4 

0.95719 

WXCL-6 

3.81003 

WXCL-7 

0.96541 

Annual 

All  Hours 

WXCL-1 

0.76201 

WXCL-2 

0.78650 

WXCL-3 

0.57538 

WXCL-4 

0.86317 

WXCL-6 

3.55605 

WXCL-7 

0.80706 

Winter 

All  Hours 

WXCL-1 

0.59869 

WXCL-2 

0.78476 

WXCL-3 

1.05942 

WXCL-4 

4.25872 

WXCL-5 

2.37126 

WXCL-6 

3.48083 

WXCL-7 

3.15858 

Annual 

All  Hours 

WXCL-1 

0.64867 

WXCL-2 

0.73165 

WXCL-3 

0.89949 

WXCL-4 

3.24614 

WXCL-5 

2.02381 

WXCL-6 

3.48151 

WXCL-7 

2.92015 

Summer 

Night 

WXCL-1 

0.88549 

WXCL-8 

0.95012 

WXCL-9 

0.74849 

Sununer 

Day 

WXCL-1 

0.85728 

WXCL-8 

1.78054 

WXCL-9 

0.67480 

Annual 

Night 

WXCL-1 

0.81240 

WXCL-8 

0.66669 

WXCL-9 

0.72328 

Annual 

Day 

WXCL-1 

0.91409 

WXCL-8 

1.22850 

WXCL-9 

0.71903 
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TABLE  9.  THE  NAMES  OF  THE  BEST  ANALYTICAL  DISTRIBUTIONS  TO 
REPRESENT  THE  OBSERVED  WIND  SPEED  DISTRIBUTIONS 


_ Time  Period  Wx  Class  Distribution 


European  Lowlands 

Winter 

All  Hours 

WXCL-1 

Gamma 

WXCL-2 

Gamma 

WXCL-3 

Weibull 

WXCL-4 

Log-Normal 

WXCL-6 

Weibull 

WXCL-7 

Gamma 

European  Lowlands 

Annual 

All  Hours 

WXCL-1 

Weibull 

WXCL-2 

Weibull 

WXCL-3 

Weibull 

WXCL-4 

Log-Normal 

WXCL-6 

Weibull 

WXCL-7 

Log-Normal 

European  Highlands 

Winter 

All  Hours 

WXCL-1 

Log-Normal 

WXCL-2 

Log-Normal 

WXCL-3 

Log-Normal 

WXCL-4 

Weibull 

WXCL-5 

Gamma 

WXCL-6 

Gamma 

WXCL-7 

Log-Normal 

European  Highlands 

Annual 

All  Hours 

WXCL-1 

Log-Normal 

WXCL-2 

Log-Normal 

WXCL-3 

Log-Normal 

WXCL-4 

Weibull 

iCL-5 

Log-Normal 

WXCL-6 

Log-Nonnal 

WXCL-7 

Log-Normal 

Mideast  Desert 

Summer 

Night 

WXCL-1 

Weibull 

WXCL-8 

Weibull 

WXCL-9 

Weibull 

Mideast  Desert 

Summer 

Day 

WXCL-1 

Weibull 

WXCL-8 

Weibull 

WXCL-9 

Weibull 

Mideast  Desert 

Annual 

Night 

WXCL-1 

Gamma 

WXCL-8 

Weibull 

WXCL-9 

Gamma 

Mideast  Desert 

Annual 

Day 

WXCL-1 

Weibull 

WXCL-8 

Weibull 

WXCL-9 

Weibull 
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EUROPEAN  LOWLANDS 
WINTER  ALL  (LST) 

WEATHER  CLASS:  WXCL-1 

GAMMA  DISTRIBUTION  (REGRESSION) 

SHAPE:  2.B0913 

SCALE:  1.42770 

NON-CALM  OCCURRENCES:  54211 

OBSERVED  DATA:  . . 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 
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Regression  fitted  gamma  distribution  for  European  Lowlands  winter 
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EUROPEAN  LOMLANDS 
ANNUAL  ALL  (LST) 

WEATHER  CLASS:  WXCL-1 

WEIBULL  DISTRIBUTION  (REGRESSION) 

SHAPE:  1. 09619 

SCALE:  4.20870 

NON-CALM  OCCURRENCES:  222626 

OBSERVED  DATA:  . 

FITTED  DISTRIBUTION:  - 
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Figure  2.  Regression  fitted  Weibull  distribution  for  European  Lowlands  wind 
speeds  for  weather  class  1. 
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I'igui'o  3.  Regression  fitted  gamma  distribution  for  European  Lowlands  winter 
wind  speeds  for  weather  class  2. 
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Figure  4.  Regression  fitted  Weibull  distribution  for  European  Lowlands  wind 
speeds  for  weather  class  2. 
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Figure  6.  Regression  fitted  Uelbull  distribution  for  European  Lowlands  wind 
speeds  for  weather  class  3. 
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Figure  7.  Regression  fitted  log-nomal  distribution  for  European  Lowlands 
winter  wind  speeds  for  weather  class  4. 
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EUROPEAN  LOWLANDS 
ANNUAL  ALL  (LST) 

WEATHER  CLASS:  WXCL-4 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  1.23490 

MEAN:  1 . 84928 

NON-CALM  OCCURRENCES:  5456 

OBSERVED  DATA: 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 
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Figure  8.  Regression  fitted  log-normal  distribution  for  European  Lowlands 
wind  speeds  for  weather  class  4. 
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EUROPEAN  LOWLANDS 
WINTER  ALL  (LST) 

WEATHER  CLASS:  WXCL-6 

WEIBULL  DISTRIBUTION  (REGRESSION) 

SHAPE:  2.63491 

SCALE:  2.32S25 

NON-CALM  OCCURRENCES:  10775 

OBSERVED  DATA:  - 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 
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Figure  9.  Regression  fitted  Weibull  distribution  for  European  Lowlands 
winter  wind  speeds  for  weather  class  6. 
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Figure  10.  Regression  fitted  Welbull  distribution  for  European  Lowlands  wind 
speeds  for  weather  class  6. 
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Figure  11.  Regression  fitted  ganuna  distribution  for  European  Lowlands  winter 
wind  speeds  for  weather  class  7. 
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EUROPEAN  LOWLANDS 
ANNUAL  ALL  (LST) 

WEATHER  CLASS:  WXCL-7 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  1.60390 

MEAN:  2.61222 

NON-CALM  OCCURRENCES:  27736 

OBSERVED  DATA:  . 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 
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Figure  12.  Regression  fitted  log-normal  distribution  for  European  Lowlands 
wind  speeds  for  weather  class  7. 
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EUROPEAN  HIGHLANDS 
WINTER  ALL  CLST) 

HEATHER  CLASS:  WXCL-1 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  a. 47344 

MEAN:  3.42020 
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Figure  13.  Regression  fitted  log-normal  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  1. 
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EUROPEAN  HIGHLANDS 
ANNUAL  ALL  (LST) 

WEATHER  CLASS:  WXCL-1 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  2.15526 

MEAN:  3.20655 

NON-CALM  OCCURRENCES:  137153 

OBSERVED  DATA: 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 
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Figure  14.  Regression  fitted  log-normal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  1. 
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Figure  15.  Regression  fitted  log-norraal  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  2. 
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Figure  16.  Regression  fitted  log-nonoal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  2. 
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EUROPEAN  HIGHLANDS 
WINTER  ALL  (LST) 

WEATHER  CLASS:  WXCL-3 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  4.20346 

MEAN:  5.69531 

NON-CALM  OCCURRENCES:  62402 

OBSERVED  DATA: 

FITTED  DISTRIBUTION:  - 

OBSERVED  -  FITTED:  . 


Figure  17.  Regression  fitted  log-normal  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  3. 
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Figure  18.  Regression  fitted  log-normal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  3. 
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WEATHER  CLASS:  WXCL-4 

WEIBULL  DISTRIBUTION  (REGRESSION) 

SHAPE:  .77226 
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Figure  19.  Regression  fitted  Welbull  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  4. 
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Figure  20.  Regression  fitted  Weibull  distribution  for  European  Highlands 
wind  speeds  for  weather  class  4. 
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Regression  fitted  gamma  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  5 . 
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Figure  22.  Regression  fitted  log-normal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  5. 
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Figure  24.  Regression  fitted  log-normal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  6. 
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Figure  25.  Regression  fitted  log-normal  distribution  for  European  Highlands 
winter  wind  speeds  for  weather  class  7. 


50 


PROB  (%)  CUM  PROB  (%) 


EUROPEAN  HIGHLANDS 
ANNUAL  ALL  (LST) 

WEATHER  CLASS:  WXCL-7 

LOG-NORMAL  DISTRIBUTION  (REGRESSION) 

SIGMA:  11.24772 

MEAN:  6.13322 
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Figure  26.  Regression  fitted  log-normal  distribution  for  European  Highlands 
wind  speeds  for  weather  class  7. 
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MIDEAST  DESERT 
SUMMER  NIGHT  (LST) 

WEATHER  CLASS:  WXCL-1 

WEIBULL  DISTRIBUTION  (REGRESSION) 

SHAPE:  1.79902 

SCALE:  4.76244 
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Figure  27,  Regression  fitted  Weibull  distribution  for  Hideast  Desert  night¬ 
time  summer  wind  speeds  for  weather  class  1. 
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Figure  28.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  day¬ 
time  summer  wind  speeds  for  weather  class  1. 
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Figure  29.  Regression  fitted  gamma  distribution  for  Mideast  Desert  nighttime 
wind  speeds  for  weather  class  1. 
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Figure  30.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  daytime 
wind  speeds  for  weather  class  1 . 
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Figure  31.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  night 
time  summer  wind  speeds  for  weather  class  8. 
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Figure  32.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  daytime 
summer  wind  speeds  for  weather  class  8. 
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Figure  33.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  night 
time  wind  speeds  for  weather  class  8. 
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Figure  34.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  daytime 
wind  speeds  for  weather  class  8. 
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Figure  35.  Regression  fitted  Weibull  distribution  for  Hideast  Desert  night 
time  summer  wind  speeds  for  weather  class  9. 
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Figure  36.  Regression  fitted  Welbull  distribution  for  Mideast  Desert  daytime 
summer  wind  speeds  for  weather  class  9. 
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rigure  37.  Regression  fitted  gamma  distribution  for  Mideast  Desert  nighttime 
wind  speeds  for  weather  class  9. 
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Figure  38.  Regression  fitted  Weibull  distribution  for  Mideast  Desert  daytime 
wind  speeds  for  weather  class  9. 
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